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Abstract. The purpose of this work is to generalize part of the theory behind Faugere's 
"F5" algorithm. This is one of the fastest known algorithms to compute a Grobner basis of 
a polynomial ideal I generated by polynomials f±, . . . , f m . A major reason for this is what 
Faugere called the algorithm's "new" criterion, and we call "the F5 criterion"; it provides 
a sufficient condition for a set of polynomials G to be a Grobner basis. However, the F5 
algorithm is difficult to grasp, and there are unresolved questions regarding its termination. 

This paper introduces some new concepts that place the criterion in a more general 
setting: 5-Grobner bases and primitive 5-irreducible polynomials. We use these to propose 
a new, simple algorithm based on a revised F5 criterion. The new concepts also enable us 
to remove various restrictions, such as proving termination without the requirement that 
fi j • • • > fm be a regular sequence. 



1. Introduction 

Since their introduction in |3j by B. Buchberger, Grobner bases and their computation 
have attracted significant attention in the computer algebra community. The best-known 
algorithm used to compute a Grobner basis is the original algorithm due to Buchberger |3j, 
and named after him. Its efficiency has been constantly enhanced through the years, but 
there remains room for improvement. Various criteria have since been introduced to detect 
useless computations - for example, (3j HJ E] - - but even so, the algorithm spends most of 
its time reducing polynomials to zero ("zero reductions"). 

In Lazard pointed out that one can view the computation of a Grobner basis as the 
reduction to row-echelon form of the Macaulay matrix of the ideal. This led to the Staggered 
Linear Basis algorithm of Gebauer and Moller in [8j , as well as the "F4" algorithm of Faugere 
in [6]. Moller, Mora, and Traverso exploited the relationship between zero reductions and 
syzygies [12], but although the algorithm they presented successfully detected many zero 
reductions, in practice it took too much memory and time (see Section 8 of [12] ). In [5], 
Faugere combined aspects of these approaches into algorithm "F5", which for a certain class 
of polynomial system eliminates all zero reductions. This algorithm exhibits impressive 
performance. 

By Faugere's admission, the theory behind the algorithm's new criterion, which we call 
the F5 criterion, is merely sketched, so as to leave more room for examples and an accurate 
description of the algorithm. The proof of the algorithm's termination and correctness were 
likewise only outlined. Additionally, some arguments were made under strong assumptions, 
such as that the input sequence fi, . . . , f m had only principal syzygies (such a sequence is 
called a regular sequence). 

We pause a moment to consider some variants of F5. Bardet described an implementa- 
tion of F5 in matrix form, where termination is ensured by manually supplying a maximal 
degree [2J. Stegers filled in some details of Faugere's proof in [TB], but stopped at two 
conjectures, one of which Gash later showed to be false |7J. 

l 
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The purpose of this paper is to present a simpler algorithm that illustrates the fundamental 
principles of F5 without sacrificing termination. We begin by defining a function S which is 
equivalent to that of Faugere, then develop a structured theory, introducing new concepts 
such as primitive S -irreducible polynomials and S -Grobner bases. These make the study of 
the problem more accessible, and suggest a new version of the F5 criterion which depends 
neither on the regularity of the input, nor on a particular ordering on the module of syzygies. 

From this theory, we develop a new, simpler algorithm. We must emphasize that the 
algorithm is a simple demonstration of the criterion, and not a deep treatment of how to 
implement a highly efficient algorithm; nevertheless, the new concepts allow us to prove cor- 
rectness and termination for any input. Note that although some F5-style algorithms provide 
explicit termination mechanisms [3 [7], these mechanisms rely on previously-developed, non- 
F5 criteria to compute explicitly a maximal degree; by contrast, the termination criterion 
used here is precisely the generalized F5 criterion used to detect useless computations. Later, 
we show that if we know that the input is a regular sequence and we use a specific ordering 
on Syz P, we can avoid all the reductions to zero. We compare the results to both F5 and 
the Staggered Linear Basis algorithm, showing how this new algorithm differs from each. 

The paper's structure is as follows. Sections I2HU cover background material; although 
most of this is relatively straightforward, an important and novel contribution of the paper 
appears at the end of Section H] with Proposition [TU The proof of that theorem leads to 
the concept of primitive iS-irreducible polynomials, from which we obtain in Section [5] a 
new characterization theorem for a Grobner basis (Theorem IT8]) . In Section [61 we use this 
characterization to formulate the new algorithm, and we prove that it terminates correctly. 
Section [7J compares this algorithm to the Staggered Linear Basis algorithm and F5, illus- 
trating the differences concretely. Section [H] describes some conclusions and possible future 
directions. 

2. Preliminaries 

Let P = k[xi, . . . , x n ] be the polynomial ring over the field k with n indeterminates, 
let fi be any admissible ordering on T", the monoid of power products over xi,...,x n : 

T n = {nr=i< i kgn}. 

Let P m be the free P-module generated by {ei,...,e m } and let fi' be any admissiblaj 
ordering on T^, the set of module terms of P m : = {tei \ t G T, / G {1, . . . , m}}. 

Fix P = (fi, . . . , f m ) G P m and let J C P be the ideal generated by P, and define 
v : P m — > I as the P-module homomorphism such that v (e«) = fi, and let Syzp = kerf, so 
that Syzp is the module of syzygies of P, LT(SyzP) C is set of leading module terms 
of Syzp, and NS(SyzP) = \ LT(SyzP) is the normal set of the syzygies of P. 

Clearly v is surjective; therefore, as a P-module, pm /syz.F ~ I. Let ip : I — > ^/syzJ 7 be the 
P-module isomorphism between them. We use the notation LT(-) for both the leading term 
of a polynomial in P with respect to fi, and the module leading term of a module element 
in P m with respect to fi'. We will use LC(/), where / is a nonzero polynomial belonging to 
/, to denote the coefficient of LT(/). 

We are interested in finding a set of polynomials G such that G is a Grobner basis for / 
with respect to the ordering /j on T". 



On what we mean by an "admissible" ordering, see the appendix. 
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Definition 1. Let 

S : I\{0} -> NS(SyzJ^) 
/ ^ LT(V>(/)), 

where LT(^(/)) is the module leading term of the normal form of i^(f) with respect to the 
ordering // on P m . 

The key idea of Faugere is to keep track of the value of S(f) for any polynomial / we will 
work with. It is however clear from the definition that the explicit calculation of S requires, 
at least, to know a Grobner basis of Syz T which is computationally expensive to compute, 
more than a Grobner basis of / itself. In fact, we will obtain S from the fact that S(fi) = e« 
(unless ti G LT(Syz J 7 )) and from other properties of S. 

3. Properties of S 

Lemma 2 (Properties of S). Let f, fi, f 2 G I \ {0}. The following hold: 

(1) If S{h) > S(f 2 ) then: 

S(f 1 + f 2 )=S(f l ). 

(2) If <S(/i) = S(f 2 ) — a an d there is no A G k* = k \ {0} such that f\ = Xf'2 ; then there 
exist a and (3 in k* such that: 

S(aft + (3f 2 ) < a. 

(3) Let t G T\ then 

S{tf)=tS{f) <=► tS(f)e NS(Syz^), 
S(tf) < tS(f) <=► tS(f) G LT(Syz^). 

Proof. ([T]) and (j2J) are trivial. 

In order to prove ([3]), let S(f ) = a = re». By the definition of S we have: 

f = v {arci + smaller terms) , 

where a G fc*, r G TP, and the argument of t> is in its normal form with respect to Syz^ 7 . 
Multiplying both sides by t, we get: 

tf = v {atrci + smaller terms) . 

If tTCi = ta G NS(Syz J 7 ), the leading term of the normal form of atrci + • • • is ta and, in 
this case, S(tf) = ta = tS(f). Otherwise, ta G" NS(Syz J 7 ), so the normal form has a leading 
term which is strictly smaller than ta and we have S(tf) < tS(f). □ 

Corollary 3. To decide whether S(tf) = tS(f), it suffices to know NS(SyzJ r ) or, equiva- 
lently, LT(SyzJ r ). Also, iftS(f) = S(g) for some g G 7\{0}, then since S(g) G NS(SyzJ r ) ; 
we can conclude that tS(f) = S(tf). 

One of the key concepts of the classic theory of Grobner bases is the polynomial reduction: 
one says that / G P \ {0} reduces with a h G P \ {0}, if there exist a G k* and t G P such 
that LT(/ - ath) < LT(/), denoted 

„ h 

where g = f — ath. 

We now introduce a special kind of reduction for a polynomial /, which takes in consid- 
eration the value of <S(/). 
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Definition 4 (S-reduction). Let f,hEl\ {0}, g E I and a E T^. We say that / S-reduces 
with respect to a to g with h, 



if there are t E T n and a E k* such that: 

• LT(g) < LT(/) and / — ath = g, and 

• S(th) < a. 

When we omit to specify a, we assume a = S(f). 

Note that this reduction is defined only for polynomials / which belong to the ideal I, 
and not for abitrary elements of the ring P. Also, when a = S(f), since S(th) < S(f) we 
have S(g) = S(f — ath) = S(f). Hence, when performing one, or more, 5-reduction steps 
with a polynomial: 



we have S(f) = S(fi) = ... = S(f k ) and LT(/) > LT(/i) > . . . > LT(/ fc ); that is, the value 



of S is kept constant, while the leading term decreases. 

Let us consider how to characterize those elements which cannot be further iS-reduced 
with respect to a given a E T^. The following definition is natural: 

Definition 5 (5-irreducible polynomial). We say that / E I is S -irreducible with respect 
to a E T™ if / = or if there is no h E I which S- reduces / with respect to a. As before, if 
we do not specify a, we assume a = S(f). Note that this definition depends on the values 
of J, T and \J . 

We could look for a criterion which decides whether a given set of nonzero polynomials 
G is a Grobner basis by looking at the values of S(g) for all g in G. However, it is wiser 
to characterize a set of polynomials with a property similar to that of a Grobner basis, but 
which also accounts for S. We therefore introduce the following: 

Definition 6 (5-Grobner basis). We say that G C I is an S- Grobner basis if for each <S- 
irreducible polynomial / E I \ {0}, there exist g E G and t E T n such that LT(ig) = LT(/) 
and S(tg) = S(f). 

Remark 7. An iS-Grobner basis depends on: 

• the ideal /, 

• the term ordering /j on T", 

• the m-tuple of generators J 7 , 

• the ordering // on T^. 

We will prove in the following section that an iS-Grobner basis is a Grobner basis in the 
usual sense. While Definition [6] is not especially useful from a computational point of view, 
inasmuch as it is quantified over an infinite set, Theorem [18] will provide us an equivalent 
criterion that is quantified over a finite set. Before we can prove it, however, we need to 
consider some properties of 5-reductions. 





4. Properties of S-reductions 



In this section we will prove the main facts which will lead to the characterization we are 
looking for. 
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Definition 8. Let 

<p: LT(J) -> NS(SyzJ-) 

t H- min{5(/) | / G /, LT(/) = t}. 

In other words, if t belongs to LT(J), ip is the minimum value S can take on a polynomial 
whose leading term is t. It follows that, for any f £ I\ {0}, </?(LT(/)) < «S(/) always holds. 

Lemma 9. (p is a bijection, and the inverse function of (p has an explicit formula: (/? _1 (cr) = 
min{tf G T n | 3/ e /, LT(/) = f , S(f) = a}. 

Proof. We show that (p is both injective and surjective. 

Injective:: By way of contradiction, suppose there exist a G NS(SyzJ r ) and ti,t 2 G 
LT(J) such that t\ > t 2 and a = tp(ti) = <p{t 2 ). Then we can find /i, f 2 G / such that 
LT(/i) = h, LT(/ 2 ) = *2, and 5(A) = «5(/ 2 ) = o". By LemmaEl there exist a,p e k* 
such that S(af\ + (3f 2 ) < o~, but LT(a/i + f3f 2 ) = t±, and therefore <p{t\) < cr, 
contradicting the hypothesis. 

Surjective:: Let cr = re^ G NS(Syz J 7 ), define 

t = mm{t' G T | 3/ G /, LT(/) = t', = a}. 

(This set is not empty because it contains LT(r/i).) Let / G / be a polynomial with 
LT(/) = t and S(f) = cr; obviously (p(t) < a. By way of contradiction, suppose that 
(p(t) < a. Then there exists /' G / such that LT(/') = t and «S(/') < cr. We can now 
choose a, a' G A;* with LT(a/ + a' f) < t such that S(af + a'f) = a. The existence 
of af + a'/' contradicts the minimality of t; therefore, (p(t) = a. 

□ 

The fact that (p is a bijection will play a crucial role in most of the subsequent proofs. 

Theorem 10 (^-reduction theorem). Let f G I and a G such that f is S -irreducible 
with respect to a and f is of the form f = v(aa + smaller terms), for some a G k* . 
Either the following equivalent propositions hold: 

(a) / = 0, 

(b) a G LT^yzJ 7 ), 

or the following equivalent propositions hold: 

(1) f*o, 

(2) a G NS^yzJ 7 ), 

(3) f^O anda = S(f)=<p(LT(f)). 

Proof. 

[3] =>■ [bl Suppose / = 0, then = / = v(aa + ••■). It follows that cr G LT(Syz J 7 ). 

[b] =^ [3} Assume by way of contradiction that a G LT(Syz^-") and / ^ 0. Let t = LT(/), 
and consider a' = <p(t) G NS(Syz J 7 ). There exists g G / such that LT(g) = t and 
5(g) = cr'; since cr' < cr, g is an (S-reductor for / with respect to cr, contradicting the 
fact that / is iS-irreducible. Therefore, / = 0. 

[TJ =>■ [2} Assume by way of contradiction that / ^ and cr G" NS (Syz J 7 ). Then a G 
LT(SyzJ 7 ). Let £ = LT(/), and consider cr' = ip (t) G NS(SyzJ 7 ). There exists 
g G / such that LT (g) = t and 5 (g) = cr'; since cr' < cr, g is an iS-reductor for / 
with respect to a, constradicting the hypothesis that / is iS-irreducible. Therefore, 
cr G NS (Syz J 7 ). 
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[2] =>- [3} Assume a G NS(SyzJ r ). Necessarily, a = S(f). Suppose now that a ^ 
</?(LT(/)); then S(f) > y>(LT(/)). Therefore there exists a polynomial g E I such 
that i = LT(#) = LT(/) and <p(t) = S(g) = <p(LT(f)) < S(f). It follows that g is an 
5-reductor of /, and / is not 5-irreducible. 

GS^CD Obvious. 

□ 

Theorem [TU] implies that it only makes sense to consider those polynomials / that are 
iS-irreducible with respect to S(f). Also, an 5-reduction yields if and only if performed 
with respect to a o G LT(Syz J 7 ); conversely, if an iS-reduction yields a non-zero polynomial, 
then we know that it was performed with respect to some a G NS(Syz^ r ). 

Remark 11. Observe that a polynomial / is 5-irreducible iff «S (/) = </?(LT(/)); oth- 
erwise, S (f) > </? (LT (/)), and we could find g G I such that LT (g) = LT (/) and 
S (g) = ip (LT (/)), so that g would 5-reduce /. 

In strict analogy with the classic Grobner basis theory we have the following result: 

Proposition 12. If G is an S-Grobner basis then for any nonzero f G I such that f is not 
S -irreducible, there exists g G G and t G T n such that: 

• LT(tg) = LT(/), 

• S(tg)=tS(g)<S(f). 

That is, it is always possible to find an S-reductor for f in G. 

Proof. Since / is not 5-irreducible, take h 5-irreducible such that LT (/) = LT (h). From the 
remark above, S (h) < S (/), so h is an 5-reductor of /. We can then find t G T and g G G 
such that t LT (g) = LT (h) = LT (/) and (using Corollary EJ S (tg) = tS (g) =S(h). □ 

This fact combined with lemma [9] leads immediately to: 

Proposition 13. If G is an S-Grobner basis, then G is a Grobner basis with respect to the 
ordering /i onP. 

Proof. For any t G LT(I), Lemma [9] implies that there exists o G such that (p~ l (o~) = t. 
Let / G / such that LT(/) = t and S(f) = a. From Proposition [T2| we may assume that / 
is 5-irreducible (if not, 5-reduce it). Then 3g G G, u G T n such that LT(ug) = LT(/) = t. 
Hence the set {LT(g) | g G G} generates LT(J) and (G) C J. Therefore G is a Grobner basis 
for /. □ 

Proposition 14. Every S-Grobner basis contains a finite S-Grobner basis. 

Proof. Let G = {gi}iei be an 5-Grobner basis, and 

i?: G P®P m 

g { ^ (LT(g i ),S(g i )). 

The image $(G) generates a P-submodule M of the P-module P © (P) m = (P) m+1 . This is 
a noetherian module; therefore, there exists a finite subset J of X such that ■&{G' S ) generates 
M, for some G' = {gj}jej- 

We claim that G' is itself an iS-Grobner basis. To see this, let / G / be an 5-irreducible 
polynomial. By definition, S (f) G NS (Syz J 7 ). Since G is an iS-Grobner basis, we can find 
a gt G G and a t G T n such that t«S = s{t gi ) = S(f) and t LT (g t ) = LT(t^) = LT(/) 
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(using Lemma [2]j3]) for tS (gi) = S (tgi)). If i G J , then G G' and we're fine. Otherwise, 
i G X \ J; since $(gi) G M, there exist j% G J7" and £j G T n such that 

This implies that = t i S{g ji ) = Sfag^) and LT(^) = LT^^-.). Then <S(%) = 

5 (£ (tiPjJ) and LT(t(yfj) = LT (i {Ug^)), so we have found gr^ G G" and £ • £j G T n which 
satisfy the iS-Grobner basis property for /. □ 

The elements of G' from the proof above will prove critically important when we examine 
our algorithm, so we will identify them by a special term. 

Definition 15 (Primitive S"-irreducible polynomial). We say that a nonzero polynomial / 
iS-irreducible with respect to S{f) is primitive S -irreducible if there are no polynomials 
/' G J \ {0} and terms t G T n \{l} such that /' is 5-irreducible, LT(t/') = LT(/) and 
S{tf') = S{f). 

The proof of Proposition [TH implies that if we have an iS-Grobner basis G, then we can 
obtain a finite 5-Grobner basis by keeping a subset of primitive 5-irreducible polynomials 
with different leading terms. Hence there exist iS-Grobner bases which contain only primitive 
iS-irreducible polynomials. 

5. The main result 

First we adapt the definition of a normal pair in [5] to reflect primitive iS-irreducible 
polynomials. 

Definition 16 (Normal Pair). Given gi,g 2 G J \ {0}, let Spol(<7i, #2) = U\g\ — u 2 g 2 be the 
S-polynomial of g 1 and g 2 ; that is, itj = ^l^^ct^) • ^ e sa y (fi^'^ 2 ) * s a norma ^ V a i r 
if: 

(1) g>j is a primitive iS-irreducible polynomial for i = 1,2, 

(2) S^ft) = LT( Ui )S( gi ) for z = 1,2, 

(3) S(u igi ) ^ S{u 2 g 2 ). 

Remark 17. With this definition, if {gi,g2) is a normal pair, then 

S{$yo\(g 1 ,g 2 )) = max{S{u 1 g 1 ),S{u 2 g2)) 

will always hold. In addition, if S{uig\) > S{u 2 g 2 ), then U\ 7^ 1, as if U\ were 1, g 2 would 
be an 5-reductor of Q\. Therefore <S(Spol(#i, g 2 )) > m8x(S{9i), £{92))- 

Theorem 18 (F5 criterion). Suppose that G is a set of S -irreducible polynomials of I , such 
that: 

• for each i = 1, . . . ,m such that ^ LT(Syz J 7 ) there exists gi G G such that S{gi) = 
Ci, and 

• for any g±,g 2 G G such that {91,92) is a normal pair, there exist g G G and t G T n 
such that tg is S-irreducible and S{tg) = S{Spo\{gi, g 2 )) . 

Then G is a S-Grobner basis of I . 

Remark 19 (Rewritable criterion). Note that the second condition does not explicitly involve 
the S-polynomial of a pair {gi, g 2 ), but cares only about <S(Spol(<7i, g 2 ))- Hence, we can think 
of this as a criterion to choose elements of NS(Syz J 7 ) instead of polynomials. Additionally, 
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if two or more normal pairs are such that S takes the same value on their S-polynomials, we 
can freely consider just one of them. 

Proof. As noted at the end of the previous section, we may, without loss of generality, assume 
that the elements of G are primitive 5-irreducible and have distinct leading terms. By way 
of contradiction, suppose that there exists a minimal a G NS(Syz J 7 ) and an iS-irreducible 
/ G I\ {0} with S (/) = a and the 5-Grobner basis property does not hold for / and a. 
That is, for all g G G and for all t G T\ LT (tg) ^ LT (/) or S (tg) ^S(f). 

The first hypothesis implies that there exist at least one primitive iS-irreducible g G G and 
some tgT™ such that rS(g) = S (/) = a; among the possible choices for g and r, pick one 
which minimizes LT(rg). By Lemma , <S (rg) = rS (g) = o. Hence LT (rg) 7^ LT(/). 
By Remark [11] S (/) = tp (LT (/)), and by LemmalU LT (rg) > LT (/). In addition, we have 
S (rg) = S (/) = ip (LT (/)) ^ cp (LT (rg)), so again by Remark ITTj rg is not iS-irreducible. 

By Lemma [2]j2]) , there exist a, (3 G k* such that S(af + j3rg) = a' for some a' < a. Since a 
was chosen to be the minimal element of NS (Syz J 7 ) such that the iS-Grobner basis property 
does not hold, Definition [6] and Proposition [12] applied to af + f3rg imply that there exist 
g' G G and r' G T" such that LT (r'g') = LT (af + f3rg) = LT (rg) and 

S (r'g') = r'S (g') < S (af + (3rg) = a' < a = S (rg) . 

Clearly g ^ g' . 

It follows that (g, g') is a normal pair. From the second hypothesis, we know that there 
exist g" G G and r" G T n such that r"g" is iS-irreducible and S(r"g") = S (Spo\(g , g')) . 
Write f Spol(<7, g') = ^yrg — 'y'r'g', for some 7, 7' G k*, where f is the gcd of r and r'. Since 
(g,g') is a normal pair and er G NS (Syz J 7 ), 

a = rS (g) = rS (Spol (<?, g')) = ?S (r" g") = S (?r"g") . 

By Remark dU S (r"g") = <p (LT (r"g")), so we have 

LT(rV) = {S(r"g")) < LT(Spol(^^'))- 
Multiplying both sides by f, we have 

LT(frV) < mirg ~ l'r'g') < LT(rg). 
The existence of g" and fr" contradicts the choice of g and r. □ 

6. The algorithm 

We shall now present a simple algorithm which computes as 5-Grobner basis of an ideal 
based on the criterion. This algorithm is quite different from Faugere's, in that it is a direct 
application of the criterion. In particular, it does not involve reductions that yield more 
then one result, nor the more rigorous simplification rules. See Section 17.21 for a detailed 
discussion. 

One first problem is that to check condition [2] of definition [TBI we need to know LT(Syz J 7 ), 
since 

S(tf) = tS(f) ^ tS(f) t LT(Syz^). 

We almost never know this before hand; therefore, we introduce a new variable L, a subset 
of LT(Syz J 7 ). At the beginning of the algorithm, we simply assume L = 0. We make use 
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of L whenever we need to check if tS(f) = S(tf) by checking whether tS(f) belongs to 
(L) C P m , the P-module generated by L. We then replace condition [2] of definition [T6l by: 

S(uigi) = LT(ui)S(gi) <^=> LT(u i )S(g i ) £ (L). 

By doing so, we end up considering more pairs than we should, but we do not skip any 
legitimate pair. 

So, when ((71,(72) is a normal pair (with the weakened condition [2]), we calculate a poly- 
nomial / = Spol((7i, (72) and a a = max(ui^i, u^g-i). Thereafter we 5-reduce / with respect 
to a. Note that a satisfies the hypothesis of Theorem [TUJ If the iS-reduction yields 0, we 
know that a G LT(Syz J 7 ); accordingly, we enlarge L by inserting a. Otherwise, we obtain a 
nonzero polynomial, which tells us that a = S(f). 

G is the set which will contain the iS-Grobner basis; we add elements to G as we find them. 
For each element g we add to G, we also store S(g); thus, G is more precisely a set of pairs 
((7,0"). When an iS-reduction returns a nonzero polynomial /, we insert / into G. Initially, 
G = 0, rather than a set containing {fi}, since we do not know if is iS-irreducible. 

B is the set of pairs of the form (/, a), where / is a polynomial that we 5-reduce with 
respect to o. Initially, we know that S(fi) = ef, therefore, we initialize B = {(/», ej)}^. 

The idea of the algorithm is to build an iS-Grobner basis by finding its elements in as- 
cending value of iS; that is, always to choose (/, a) G B such that a is minimal. (See step 

M) 

Remark 20. In virtue of remark [T9| for each a we can keep in B at most one polynomial / 
such that S (/) = a. For the same reason we can, at any time, remove (/, a) from B if we 
can find another polynomial /' such that S (/') = a and LT (/') < LT (/). 

In practice we will remove from B a pair (/, a) if we can find at 6 P and a (/', a') G GUB 
such that to' = a and t LT (/') < LT (/). 

The pseudo code of the algorithm is the following: 

Algorithm 21. 

Input: T = (/1, . . . , f m ): an element of P m , 

fi: an ordering on T n , 

fi' : an ordering on T^. 
Output: G: an S-Grdbner basis of I = (f x , . . . , f m ). 

(1) L:=0 

(2) G := 

(3) B :={(/!, ei ),...,(/ m ,e m )} 

(4 ) While B^0 

(a) B:={(f,a)eB\a?(L)} 

(b) Remove from B any (/, a) such that we can find (f',a') G G U B, t G T n 
satisfying to 1 = a and LT (tf) < LT (/) 

(c) Pick (/, a) G B with minimal a. 

(d) f := S-reduce(f,a,G) 

(e) If f j£ then 

(i) B := UpdatePairs(L, G, B, (f, a)) 

(ii) G:=GU{(f,a)} 

(f) Else 

(1) L :=LU{(i} 
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(5) Return {g : (g, a) G G} 

Note that, since L may change during each iteration, some pairs we assumed to be normal 
turn out not to be normal. We remove those in step |4al 

In step |4b] we implement the idea presented in Remark [201 Note that this is an optimiza- 
tion; the algorithm will successfully terminate without this line. 

We still have to describe the two procedures Algorithm |2"T1 invokes. The first is S-reduce: 

Algorithm 22 (iS-reduce). 

Input: /; an element of I , 
a: an element ofY™ n , 

G: a set that contains the elements (g, S (g)) of an S-Grdbner basis with S(g) < a. 
Output: /; an S -irreducible polynomial with respect to a. 

(1) /:= //LC(/) 

(2) While 3 (g, S (g)) G G,t G T n such that t LT (g) = LT (/) and tS (g) < a 

(a) f:=f-tg/LC(g)) 

(b) If f = then Return 

(c) /:=//LC(/) 

(3) Return f 

This algorithm takes as input a polynomial / and a a G and, as long as there is an 
5-reductor for / in G, performs iS-reduction steps. Because of the hypothesis on G we know 
we obtain an iS-irreducible polynomial with respect to a. 

The second is UpdatePairs: 

Algorithm 23 (UpdatePairs). 

Input: L: a subset of LT (Syz J 7 ), 

G: a set that contains the elements (g, S (g)) of a S-Grdbner basis with S (g) < a, 
B: a set that contains elements (g, a g ) of polynomials that have yet to be considered, 
(f,S(F)): where / G /\ {0}. 

Output: B' : a set of pairs (/', a) that satisfy Theorem [Td[ produced by the criterion. 

(1) B' := 

(2) For each (g,S (g)) G G, if (f,g) is a normal pair 

(a) Compute U\, u<i such that Spol (/, g) = U\f + Uig 

(b) a := max (« x 5 (/) , u 2 S (g)) 

(c) B':=B'U{(Spo\(f,g),a)} 

(3) Return B' U B 

Proposition 24. Alaorithml21\ terminates. 



Proof. First we show that step 4(f)i is executed only a finite number of times. 



Because of step HU at a given time, we only consider a that do not belong to L; so when 



we execute step 4(f)i we really enlarge the P-module generated by L. Since P m is noetherian 



this can happen only a finite number of times. 



Also, that step 4(e)i is executed only a finite number of times. First note that if / is not 
primitive iS-irreducible (that is, / is only iS-irreducible) , then Algorithm [231 does nothing, so 
no new polynomials are generated. In the proof of Proposition [TU we see that an 5-Grobner 
basis contains only a finite number of primitive iS-irreducible polynomials. This completes 
the proof. □ 
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Theorem 25. Alaorithml21\ computes an S-Grdbner basis of I. 

Proof. This is a direct consequence of criterion of Theorem ITB1 Previous remarks have shown 
that G contains only 5-irreducible polynomials, and the initial value of B ensures that the 
algorithm satisfies the first condition. For the second condition, for each normal pair (g\, g-z), 
we ensure that we have a polynomial / and a monomial t such that S(tf) = <S(Spol(<7i, g%))- 

□ 

The fact that non-primitive iS-irreducible polynomials do not generate any new pairs plays 
a central role in this proof of termination. Without it, the thesis does not hold: if we drop 
condition (JTJ) of Definition [161 it is possible that the algorithm could enter an infinite loop, 
computing an infinite number of polynomials of the form tif where {ti} is an infinite set of 
terms and / is an iS-irreducible polynomial and each of the tif is 5-irreducible itself. (This 
occurs, for example, in the implementation of [15J.) 

7. Comparison with previous work 

In this section, we consider how this algorithm is both similar and different to two al- 
gorithms in past work: the staggered linear basis algorithm of Gebauer and Moller [8J and 
the F5 algorithm of Faugere j5]. We also illustrate explicit differences on three particular 
examples. 

7.1. Comparison with Staggered Linear Bases. The Staggered Linear Basis algorithm 
(in the rest of this section, SLB) |8J introduced a special kind of Grobner basis. 

Definition 26. The set B C I is a staggered linear basis of the ideal I if for all f G P 

• if /, g G B and LT (/) = LT (g), then f — g; and 

• if t G T and tf G B, then / G B. 

A full review of SLB is beyond the scope of this paper, but it is worth comparing to the 
present algorithm because both use trivial syzygies to detect zero reductions. To facilitate 
the explanation, we temporarily adopt the notation ti = LT (fi) and tjj = 1cm (ti, tj). 

SLB tracks monomial ideals for each polynomial among the generators. Initially, we have 

Zi = (ti, t%, . . . , ti-i) . 

Critical pairs (fi, fj) (with i < j) are rejected whenever G Zi. If instead the 5-poly- 

nomial of (fi, fj) is computed, then Zj is expanded by adding the ideal generated by Uj/tj. 
If reduction of the ^-polynomial results in a new polynomial being added to the basis, 
SLB also creates a new ideal 

Zk = (Zj + (t^) : (tij/tj) + (t%, . . . , . 

Despite the use of principal syzygies in the initial definition of Zj, a fundamental difference 
between the algorithms lies in the fact that SLB does not compute, let alone consider, the 
leading module term S (f) of any polynomials. So a polynomial can be iS-irreducible even if 
it is top-reducible, and the normal pairs of the F5 Criterion are not the same as the critical 
pairs of SLB. As a result, the approach in SLB behaves quite differently, and fails to detect 
certain zero reductions detected by F5 and the present algorithm. 
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7.2. Comparison with F5. At first glance, algorithm [211 may appear very different from 
the F5 algorithm. However, if we define \J to be 



i < j or 

i — j and t < s 



for any t,s E T n and 1 < i,j < m, there is an interesting relationship between 5-Grobner 
bases and LT(Syz J 7 ). Define, for 1 < / < m, tti : P m — > P as the projection on the l-th 
component, then 

(7.1) 7r,(PSyzJ0 = (/i,...,/j-i) 

where PSyz(J r ) is the P-module of principal syzygies, defined as PSyz(J r ) = (fiCj — jjtijp C 
P m ; PSyz(J-") is clearly a P-sub module of Syz(J-"). 

Suppose we have f £ I \ {0} and we know S(f) = tei, for some t G T n . It follows from 
the definition of S that / G (/i, . . . , fi). Hence, (17. ip implies that 

LT(f)e i+j G LT(PSyz(J r )) for some j > 1. 

With this choice for //, we can improve the performance of Algorithm [5T] by adding an 



instruction right after step 4(e)ii 



L := L U {LT(/)e l+1 , LT(/)e i+2 , . . . , LT(/)e m }, 

where a = tei for some t G T n . In other words, whenever we find a new element of G, we 
also find new elements of L. 

Also, due to the ordering on P m , the structure of an iS-Grobner basis G is very special. 
We find the elements of G in ascending value of S: we first find all the elements g such 
that S(g) = tei f° r some t G T n , then those g such that S(g) = tei f° r some t G T n and 
so on. Is easy to see that the real value of f] is never considered in any computation, until 
the algorithm has finished producing all the elements of G with S(g) = te, for some t G T n 
and i < I. If we make the further assumption that Syz(J r ) = PSyz(J-"), we conclude that 
the algorithm never reduces a polynomial to 0, since we discover every leading term of the 
syzygies in advance. 

Therefore, we may say that, in this case, Algorithm |2~T1 is incremental, as it first produces 
an 5-Grobner basis of (fx), then an 5-Grobner basis of (fi, fi) and so on, and avoids all the 
reductions to zero; this behavior is the same as Faugere's F5 algorithm. 

We can couch the use of "simplification rules" in F5 [5, sect. 6], also called the rewritable 
criterion, in vocabulary similar to that used in this paper: F5's algorithm to compute S- 
polynomials (SPol) discards any (tei, f) when 

• there exists some other (t'ei, /') G G U B U B' such that t | t', and 

• /' was computed before f. 

This concept is related to Remark [TH] in this paper; roughly we know we can "decide" how 
to obtain an iS-irreducible polynomial with a given signature. We prefer to start with a 
polynomial with the smallest leading term we know of, while in F5 just the first generated 
polynomial is kept. 

This parallel carries over to the computation of L, which here is used to prevent the 
computation of any tei G LT (Syz J 7 ) more than once. When a polynomial is reduced to zero 
in F5, the simplification rule is added even though the polynomial is discarded, and this rule 
ensures that any polynomial / with S (f) = t! Ci, where t \ t! , is not computed. In other 
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Number of zero reductions 



Algorithm 


MMT92 


Cyclic-5 


Cyclic-6 


Katsura-5 


Staggered Linear Basis 


3 


46 


446 


10 


F5 








16 





Algorithm |2"T] 








8 






TABLE 1. Number of zero reductions during execution of algorithms SLB, F5, 
and Algorithm I2T1 



Size of basis 



Algorithm (size of red. GB) 


MMT92 


Cyclic-5 


Cyclic-6 


Katsura-5 


Staggered Linear Basis 


8 


38 


99 


22 


F5 


10 


39 


202 


30 


Algorithm 1211 


10 


39 


155 


30 



TABLE 2. Size of the Grobner basis computed by algorithms SLB, F5, and Algorithm [2T1 



words, F5 has an implicit provision for avoiding the computation of non-trivial syzygies, like 
the algorithm here. 

7.3. Concrete examples. We examine how all three examples perform on three "standard" 
systems: 

• the system "MMT92", F = {yz 3 — x 2 t 2 , xz 2 — y 2 t, x 2 y — z 2 t} from [5] (this seems first 
to appear in non-homogenized form in |12]): 

• the homogenized Cyclic-5 system; and 

• the homogenized Katsura-5 system. 

We consider 

(1) the number of zero reductions; and 

(2) the size of the Grobner basis generated. 

The tests were carried out in unoptimized implementations of each algorithm in Sage |16[ (TJ 
[HI [13], and are available online. 

The results are in Tables [T] and |5J Table[T]shows that the Staggered Linear Basis algorithm 
computes some zero reductions even though the systems are regular sequences. Neither F5 
nor Algorithm [5T] computes any zero reductions except in Cyclic-6, which is not a regular 
sequence. In that system, Algorithm [21] computes a smaller basis, and it computes fewer 
zero reductions. This appears to be due to the fact that it proceeds by ascending signature 
(line He]) rather than by ascending 1cm (compare to algorithm Spol in [5]). 

8. Conclusions and future work 

This paper has reformulated the F5 criterion, which in its original form is due to [5], 
and provided a new proof of this criterion's correctness. We have introduced the ideas of 
S -Grobner basis and S -irreducible polynomials, and have shown that if a set of polynomials 
G satisfies the F5 criterion, then G is an 5-Grobner basis and not just a Grobner basis. In 
this new setting, we were able to drop many restrictions present in [5j: we can freely choose 
any ordering on P m , and there is no need to for the sequence . . . , f m ) to be regular. 
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Our statement of the criterion is quite different from the original: we require that all the 
polynomials in the set G be iS-irreducible; we require that if (jL NS(Syz J 7 ), then there exist 
g G G such that S (g) = ef, and we impose a condition on the signature S (Spol (gi, #2)); 
rather than the usual condition that 

#G 

Spol (5-1, g 2 ) = ^2 h i9i sucn that hi ^ LT LT {gi) < LT (Spol (5-1, 5-2)) • 
i=i 

(Faugere calls this latter condition (Spol {g±, g 2 ))-) We also changed the definition of normal 
pair by adding a new condition: the fact that we can consider only primitive 5-irreducible 
polynomials. 

We then proposed a simple algorithm to show an application of the new criterion. The 
algorithm presented here is mainly demonstrative, and does not include many "obvious" 
optimizations such as holding off on the computation of a new polynomial / until it is 
actually needed in step |4d] of Algorithm [2TJ 

The authors would like to thank the referees for helpful and instructive comments that 
improved the paper. 

Appendix 

In October 2011, Vasily Galkin of Moscow State University contacted us with a question 
about Proposition [HJ His question was sparked by the definition of // as an admissible 
ordering on T^. Apparently, we used the wrong word; our reference for the notation (which, 
it amazes us to report now, we did not include in the bibliography) was [10J. (As far as 
we can tell, it is the only textbook that uses the word "monomodule".) This text does not 
define an "admissible" ordering for a module; it defines either a module ordering (p. 54) 
or a compatible ordering (p. 55). It seems that when we wrote "admissible", we meant 
"compatible". Indeed, if the ordering is not compatible, the iS-Grobner basis may be infinite, 
as the following example shows. 

Example 27. Let < be the degrevlex ordering with x > y > z. Let <' be the module 
ordering with 

xe\ <' yei <' xe 2 <' ye 2 <' ze\ <' ze 2 
extended to all other module terms in the following way: tej <' uej if 

• degt < degu (total degree), or 

• deg t = deg u and 

- deg 2 1 < deg 2 u, or 

— deg 2 1 = deg z u and 

* i = 1 and j = 2, or 

* i = j and deg^ t < deg y u, or 

* i = j, deg y t = deg y u, and deg x t < deg x u. 

It is routine to verify that this is a module ordering; it is obviously not compatible with 

<. 

Let 

fi=x 2 + xy, f 2 = X y + z 2 . 

We have 



(8.1) LT (Syz ({/1, / 2 })) = {t ■ LT {f 2 e, - he 2 ) :teT n } = {t- z 2 e, :teT n } 
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For % > 3, let f { be the S'-polynomial of fi and fi_ x . We have 

r 2 2 

J3 = xy - XZ 
U = xy 3 + x 2 z 2 

f i = xy i - 1 + (-l) i x i - 2 * 1 . 

For % > 2, S (fi) = x i_2 e 2 = min {te 2 : t G T n , deg t = i — 2}. So /; is not 5-reducible by / 2 , 
. . . , even though it is top-reducible by them. Additionally, LT (/i) f LT (/,), so /i is 
iS-irreducible. 

It remains to see if the fi are primitive. Let g G I\ {0} such that LT (tg) = LT (/j) 
for some t G T n \{l}. The factors of LT (fi) are x and y on/y, so t — x a y b for some 
a, b G N. We also want S (tg) = S (fi); we claim that this is possible only if t = x a for 
some a G N. To see why, assume that S (tg) = S (/). If t«S (g) ^ S (g), we must have 
tS (g) G LT (Syz ({fi, / 2 })) ; which would imply that tS (g) = uz 2 e 1 for some u G T n . Since 
S (g) G NS (Syz ({/i, /2})), we infer that 2 > deg 2 5 (g) = deg 2 (t«S (g)), a contradiction to 
tS (g) = uzW Thus, *5 ((?) G NS (Syz ({/i, / 2 })), and tS (g) = S (tg). So tS (g) = S (tg) = 
S (fi) = x l ~ 2 e 2 . As claimed, t = x a for some a G N. Since t ^ 1, a > 0; since t | LT(/j), 
a < 2. Hence, t = x and LT (g) = y' 1 ^ 1 , but a computation of the Grobner basis shows that 
y 1 " 1 <£ LT(7) = (t/ 2 z 2 ,xz 2 ,x 2 ,xy). Hence, there do not exist g G I\ {0} and i G T n \{l} 
such that LT (^) = LT (/i) and S (tg) = S (fi). 

It follows that any iS-Grobner of I with respect to < and <' must have all the polynomials 
f 2 , fs, ■ ■ ., and is thus infinite. 
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